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The present generation of weak lensing surveys will be superseded by surveys run from 
S ' space with much better sky coverage and high level of signal to noise ratio, such as 

SNAP. However, removal of any systematics or noise will remain a major cause of 
concern for any weak lensing survey. One of the best ways of spotting any undetected 
source of systematic noise is to compare surveys which probe the same part of the sky. 
In this paper we study various measures which are useful in cross correlating weak 
lensing surveys with diverse survey strategies. Using two different statistics - the shear 
components and the aperture mass - we construct a class of estimators which encode 
such cross-correlations. These techniques will also be useful in studies where the entire 
^^ ■ source population from a specific survey can be divided into various redshift bins to 

study cross correlations among them. We perform a detailed study of the angular size 
dependence and redshift dependence of these observables and of their sensitivity to 
the background cosmology. We find that one-point and two-point statistics provide 
complementary tools which allow one to constrain cosmological parameters and to 
obtain a simple estimate of the noise of the survey. 
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.. ', 1 INTRODUCTION 

>. 
•i-H . Detection of weak lensing signals by observational teams (e.g., Bacon, Refregier & Ellis, 2000, Hoekstra et al., 2002, Van 

^^ ' Waerbeke et al., 2000, and Van Waerbeke et al., 2002) has led to a new avenue not only in constraining the background 

JH , dynamics of the universe but in probing the nature of dark matter and dark energy as well. To do so one compares observations 

. . . i with theoretical results obtained from simulations or analytical methods. 

Numerical simulations of weak lensing typically employ ray-tracing techniques as well as line of sight integration of cosmic 
shear (e.g., Schneider & Weiss, 1988, Jarosszn'ski et al., 1990, Wambsganns, Cen & Ostriker, 1998, Van Waerbeke, Bernardeau 
& Mellier, 1999, and Jain, Seljak & White, 2000,Couchman, Barber & Thomas (1999)) and provide valuable insights. 

On the other hand, analytical techniques to study weak lensing include perturbative calculations at large angular scales 
(e.g., Villumsen, 1996, Stebbins, 1996, Bernardeau et al., 1997, Jain & Seljak, 1997, Kaiser, 1998, Van Waerbeke, Bernardeau 
& Mellier, 1999, and Schneider et al., 1998) and techniques based on the hierarchical ansatz at small angular scales (e.g., Fry 
1984, Schaeffer 1984, Bernardeau & Schaeffer 1992, Szapudi & Szalay 1993, 1997, Munshi, Melott & Coles 1999). Ingredients 
for such calculations include Peacock & Dodds (1996)'s prescription (see Peacock & Smith (2000) for a more recent fit) for the 
evolution of the power spectrum or equivalently the two-point correlation function. Recent studies have shown an excellent 
agreement between analytical results and numerical simulations of weak lensing effects (Valageas 2000a & b; Munshi & Jain 
2000 & 2001; Munshi 2000; Bernardeau & Valageas 2000; Valageas, Barber & Munshi 2004; Barber, Munshi & Valageas 2004; 
Munshi, Valageas & Barber 2004). 

However, to the weak lensing effects themselves one must add various sources of noise, such as the intrinsic ellipticity 
distribution of galaxies, shot-noise due to the discreet nature of the source galaxies and finite volume effects due to finite 
survey size. There is a need to incorporate all these in a consistent and systematic way in any realistic design of survey 
strategy. A detailed formalism to tackle such issues was developed by Munshi & Coles (2002) following Schneider et al. (1998). 
Recently Valageas, Munshi & Barber (2004) have used such techniques and have incorporated the redshift distribution of 
sources to model errors in future surveys such as SNAP. 

In the near future we expect to see many weak lensing surveys probing overlapping areas of the sky with different 
survey strategies. These ground based observations will finally be superseded by surveys run from space with much better 
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sky coverage such as SNAP. Nevertheless, removal of systematics and control of the noise from individual surveys will remain 
a major cause of concern. It is therefore of utmost importance that any undetected systematics be spotted by comparing 
different surveys which probe the same regions of the sky. Such a cross comparison will finally lead to the construction of 
large maps by combining smaller maps in an optimum way. 

In our present study we focus on cross correlations among various surveys or different redshift bins in the same survey, 
with partial or complete overlap on their sky coverage. Using two different classes of statistics such as the smoothed shear 
components and the aperture mass, we construct estimators which encode such correlations. Based on analytical results from 
Munshi & Coles (2002) we generalize a recent study by Valageas, Munshi & Barber (2004) to take into account multiple 
surveys. We display the power of such techniques using SNAP class surveys where the source population can be subdivided 
into various redshift bins to study cross correlation among subsamples. We study the angular size dependence and redshift 
dependence of these classes of statistics. We also describe the dependence on background cosmological parameters. 

This paper is organized as follows: in section 2, we describe our notations regarding weak lensing observables in general 
and various filters used to smooth the data. In section 3, we introduce weak lensing estimators which take into account various 
realistic sources of noise. Specific survey geometries based on SNAP class experiments are introduced in section 4 and we 
compute the noisy cumulant correlators of the aperture mass and of the shear components. Section 5 is devoted to a discussion 
of our results. 



2 NOTATIONS AND FORMALISM 

Let us first recall our notations. Weak-lensing effects can be expressed in terms of the convergence along the line-of-sight 
towards the direction on the sky up to the redshift z 3 of the source, k(i}, z 3 ), given by (e.g., Bernardeau et al. 1997; Kaiser 
1998): 

K tf,Z.) = ^ r 3 dxw(x,Xs)5(x,V#), with 4ri-f ^^ (1 + *), (1) 

where z corresponds to the radial distance x an d T> is the angular distance. Here and in the following we use the Born 
approximation which is well-suited to weak-lensing studies: the fluctuations of the gravitational potential are computed along 
the unperturbed trajectory of the photon (Kaiser 1992). Thus the convergence k($, z s ) is merely the projection of the local 
density contrast S along the line-of-sight. Therefore, weak lensing observations allow us to measure the projected density field 
k on the sky (note that by looking at sources located at different redshifts one may also probe the radial direction). In practice 
the sources have a broad redshift distribution which needs to be taken into account. Thus, the quantity of interest is actually: 

«($) = / dz s n(z s )K('d,z s ), with / dz s n(z s ) — 1, (2) 

where n(z 3 ) is the mean redshift distribution of the sources (e.g. galaxies) normalized to unity. From eq.(l), the convergence 
k associated with a specific survey also reads: 

k(#) = / dxw(x) 5(x,T>tf), with w(x) = —^- \ dz s n(z s ) w(x,X*), (3) 

where z max is the depth of the survey (i.e. n{z s ) = for z 3 > z max ). By working with eq.(3) we neglect the discrete effects due 
to the finite number of galaxies. They can be obtained by taking into account the discrete nature of the distribution n(z 3 ). 
This gives corrections of order 1/N to higher-order moments of weak-lensing observables, where N is the number of galaxies 
within the circular field of interest. In practice N is much larger than unity (for a circular window of radius 1 arcmin we 
expect N <; 100 for the SNAP mission) therefore in this paper we shall work with eq.(3). 

Thus, in order to take into account the redshift distribution of sources we simply need to replace 3f2 m /2w(x, Xs) m e< l-(l) 
by w(x)- Therefore, all the results of Munshi et al. (2004) remain valid. Then, usual weak-lensing observables can be written 
as the angular average of k($) with some filter U: 

X = di>Ux(#) K0). (4) 

For instance, the filters associated with the smoothed convergence k 3 , the smoothed shear 7,, and the aperture- mass M ap are 
(Munshi et al. 2004): 

__ e and t7 Map = ___ 9 M__ _ 



U..=^^, U ls =-^pie^ and C7 Map = ^^ 9 ( 1 - ^ ± - ^ ) , (5) 



where are Heaviside functions with obvious notations and a is the polar angle of the vector {>. The angular radius 6 3 gives 
the angular scale probed by these smoothed observables. Note that the smoothed shear 7 S depends on the matter located 
outside of the cone of radius 9 3 . However, in practice one directly measures the shear 7(1?) on the direction $ (from the 
ellipticity of a galaxy) rather than the convergence k and 7 S is simply the mean shear within the radius 6 3 . For M ap we shall 
use in this paper the filter (5), as in Schneider (1996), but one could also use any compensated filter with radial symmetry. 
As described in Munshi et al. (2004), the cumulants of X can be written as: 
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;*">.= / ||'i\-«-(\,) / Mch/.r.vt,/,) c„( \\ { , g^ ,..., g^ ). ((,) 



i=l ^ fcl V 



or equivalently we can write: 

poo P /■ P 

(X p )c= I T\d Xt w( Xl ) T\d^W x (k ±J V 3 e s ) |||e**n«w) (5(k 1 )...5(k p )) c . (7) 



y n d »*(»)y n^wx^-Dje,) me^j wci)...^^ 



We note (..} the average over different realizations of the density field, £ p is the real-space p— point correlation function of the 
density field £ p (xi, .., x p ) = (5(xi)..<5(x p )) c , fcii is the component of k parallel to the line-of-sight, kx is the two-dimensional 
vector formed by the components of k perpendicular to the line-of-sight and Wx(k±T>9 s ) is the Fourier transform of the 
window Ux'- 

W x (k±V9 s ) = fdtiUx&e^-™. (8) 

In particular, for the smoothed convergence k s , the smoothed shear 7 S and the aperture-mass M ap we have (Munshi et al. 
2004): 

W Ks (k ± V6 s )= kM , W 7 .(k x Pg.)= k±V6s e* and M/ Map (k ± Pg s ) = ^-^ , (9) 

where is the polar angle of kx and J„ are Bessel functions of the first kind. The real-space expression (6) is well-suited to 
models which give an analytic expression for the correlations £ p , like the minimal tree-model (Valageas 2000b; Bernardeau 
& Valageas 2002; Barber et al. 2004) while the Fourier-space expression (7) is convenient for models which give a simple 
expression for the correlations {8(ki)..8(k p )) c , like the stellar model (Valageas et al. 2004; Barber et al. 2004). In the present 
study we shall also be interested in two-point cumulants such as: 

(X{X q ) c = / d X (2^Y +q - 1 w{wl / Yl dk ±J WxCk^VOs) (5(kxi) . . . S(k ±p+q )) c , (10) 

Jo J j=1 

which describe the cross-correlations between two surveys or subsamples. Here the subscripts "1,2" refer to the two surveys 
we cross-correlate. In eq.(10) we used the fact that the correlation length is much smaller than cosmological scales and the 
Dirac factor <$d(A;||i + . . . + ku p + q ) has been factorized out of (<5(kxi) . . . S(k± p+q )) c . The two-point correlators are similar to 
projected cumulant correlators studied in previous works (Szapudi & Szalay 1993; Munshi et al. 1999 ). However, unlike those 
projected cumulant correlators these objects are defined here in the zero angular separation limit where only the projection 
along the radial direction is different for the two surveys. Note that for observables like the aperture-mass which are very 
localized (long wavelengths are damped by the compensated filter) it is expected that projected cumulant correlators will only 
have non-negligible values if the angular separation is not much larger than the angular radius of the window. 
From the one-point cumulants (6)-(7) it is convenient to define the parameters S p and the coefficients t p : 

The quantities S p apply to one subsample while t p apply to two subsamples. The coefficients S p of weak lensing observables are 
closely related to the usual parameters S p which describe the departures of the underlying 3-d density field from Gaussianity 
(e.g., Valageas et al. 2004, Munshi et al. 2004). In particular, for the case of the smoothed convergence k s we have S p 3 ~ S p 
where k s is the smoothed convergence normalized by a suitable factor and S is the 3-d density contrast (e.g., Barber et al. 
2004). The normalization of the coefficients S p is such that they have a finite limit in the quasi-linear regime and show at most 
weak variations in the highly non-linear regime (but they exhibit a sharp change at the transition). Thus, from the one-point 
quantities S p we obtain the deviations from Gaussianity associated with each survey or subsample. Together with the second- 
order moments (X 2 ) (i.e. variance) they can be used to measure cosmological parameters as well as the non-Gaussianities 
built by the non-linear gravitational dynamics. In this regard, the redshift binning of the source population can also be very 
useful as it allows one to constrain the time evolution of cosmological distances and growth factors through the ratios t p . 
Obviously, the weak lensing effects associated with different redshift bins are correlated since their lines of sight probe the 
same density fluctuations at low z (where they are largest), for a fixed angular patch on the sky. In order to measure these 

as: 

<*f*2>c - - , {X 1 X 2 ) C 

m P artlcular r n= „'i„,li,v (12) 



rp+q\p/(p+q) / yP+9\?/(P+?) 



(X p q ) p/{p q> (X ) q (X 2 ) (X 2 ) 

The quantities r pq correspond to the two-point cumulants {XfX q ) c normalized in such a way that most of the dependence on 
cosmology and gravitational dynamics cancels out. Thus, if the two subsamples are highly correlated we have r pq ~ 1 while 
r pq ~ if they are almost uncorrelated. 

One often uses cumulant correlators to study the cross-correlations as a function of the separation angle 812 between two 
different angular patches. In the numerical computations we perform in this paper we shall restrict ourselves to the case of 

© 0000 RAS, MNRAS 000, 000-000 



4 Munshi & Valageas 




Figure 1. A schematic representation of cross-correlating two different source populations represented by mean source redshifts z s i and 
z s 2- Non-overlapping source distributions are considered. Survey areas are Oi and C2 respectively, the common smoothing angular radius 
is 8 S . The contributions to the scatter of weak lensing estimators due to the galaxy intrinsic ellipticities of each sample can be described 
by the quantities pi(9 B ) and P2(S S ) defined in eq.(16). 

zero angular separation. Hence we focus on the correlations between two surveys which overlap on the sky but which have 
different redshift distributions of sources (see Fig. 1). On the other hand, note that the signal would decrease for nonzero 
angular separations. 

3 LOW-ORDER ESTIMATORS AND THEIR SCATTER 

The expressions obtained in the previous section describe weak lensing effects due to the fluctuations of the matter density 
field. However, in order to handle realistic data one needs to take into account the observational noise. A first source of scatter 
is due to the finite size of the surveys. A second source of noise is merely due to the intrinsic ellipticity of galaxies, which 
cannot be avoided. As in Munshi & Coles (2003) and Valageas, Munshi & Barber (2004), this leads us to define the estimators 
M v for low-order moments: 



M =M£ 



(il,...,ip) 



i p 6jj 



with (N) p = N(N-i)..(N-p + l) = 



N\ 



(N-p)V 



(13) 



where N is the number of galaxies in the patch of size ir8 s and p is the order of the moment. Here we used the fact that 
the aperture- mass defined from the convergence n by the compensated filter Um* p given in eq.(5) can also be written as a 
function of the tangential shear 7t as (Kaiser et al. 1994; Schneider 1996): 



M, 



d0QM„ p (0)7t(0) with Q A / ap (0) 



9(0 < 9s) c / 



n6 2 



(B 



o 2 



while for the smoothed shear components 7i S (with i = 1, 2) we simply have: 
7ls = / d0Q 7js (0) 7i(0) with Q 7ls (0) = e( ^^/ s) . 



(14) 



(15) 



Thus, in eq.(13) we wrote Qj — Qxi'&j) and we use the tangential ellipticity (for M ap ) or the i-component of the ellipticity 
(for 7i s ). Indeed, in the case of weak lensing, k<1, the observed complex ellipticity e is related to the shear 7 by: e = 7 + e*, 
where e, is the intrinsic ellipticity of the galaxy. Finally, the sum in eq.(13) runs over all ordered sets of p different galaxies 
among the N galaxies enclosed in the angular radius 8 S . 

Then, neglecting any cross-correlations between the intrinsic ellipticities of different galaxies and between the matter 
density field and the galaxy intrinsic ellipticities, the estimators M p introduced in eq.(13) are unbiased estimators of low-order 
moments of the aperture-mass or of the shear components. Moreover, their dispersion a 2 (M p ) = {M 2 } — (Af p ) 2 only involves 
the variance a 2 = (|e*| 2 ) of the galaxy intrinsic ellipticities (even if the latter are not Gaussian). The latter enters the scatter 
a 2 through the combination p defined by: 

2N(X 2 } 



o 2 Gy 



with Gx — 7r6 



d0' 



x(0) , whence G A / ap = -, G liB = 1. 



(16) 



Here N is the number of galaxies within the circular field of angular radius 6 3 . The quantity p measures the relative importance 
of the galaxy intrinsic ellipticities in the signal. They can be neglected if p 3> 1. Thus, the mean of the estimator M2 and its 
scatter are given by (with (X) = 0, see also Valageas et al. 2004): 
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(M 2 ) = (X 2 ) c and a 2 {M 2 ) = {X 4 ) c + (X 2 ) 2 C 2 



1 + 
P 



r 2 



(17) 



In order to obtain eq.(17) we have averaged i) over the intrinsic ellipticity distribution, ii) over the galaxy positions and 
iii) over the matter density field, assuming these three averaging procedures are uncorrelated. Any Gaussian white noise 
associated with the detector can be incorporated into these expressions by adding a relevant correction to a 2 (whence to p). 
Note that the scatter of low-order estimators involves higher-order cumulants (up to twice the order of the moment one wishes 
to evaluate). Here we have neglected the terms due to the finite number of galaxies which scale as 1/iV and are negligible for 
practical purposes (N ^ 1). The dispersion a 2 obtained in eq.(17) takes into account two sources of noise: the galaxy intrinsic 
ellipticities (the terms which depend on p) and the finite size of the survey (the remaining terms) . 

The estimators M p defined in eq.(13) correspond to a single circular field of angular radius 9 S containing N galaxies. In 
practice, the size of the survey is much larger than 9 S and we can average over N c cells on the sky. This yields the estimators 
X p defined by: 

N c 

M P = ^-Y]M^\ whence (M p ) = (M p ) = (X p ) and a(M p ) = I^MA witri a 2 {M p ) = {M p ) - (M p ) 2 , (18) 

n — 1 
(n) 

where M p is the estimator M p for the cell n and we assumed that these cells are sufficiently well separated so as to be 
uncorrelated. 

The estimators M p and X p provide a measure of the moments (X p ) of weak lensing observables, which can be used 
to evaluate the S p parameters defined in eq.(ll). However, as shown in Valageas et al. (2004), it is better to first consider 
cumulant-inspired estimators H p . Thus, for the aperture mass where we shall restrict ourselves to third-order cumulants we 
define: 

N C 

H 3 = M 3 - 3X 2 Mi and H 3 = j ^ H { 3 n) , whence {Hz) = <M a 3 p ) c , (19) 



u=l 



while for the shear components which are even quantities (so that all odd-order moments vanish) we need to go up to fourth 
order and we define: 



H i = M i -6M 2 M 2 + 3M 2 and Ha = -£- ^ H T' , whence (Hi) = (it s )c- (20) 

n=l 



The interest of H 3 and Hi is that their scatter is smaller than for M3 and Mi, see Valageas et al. (2004). Besides they directly 
yield the one-point cumulants (here we neglected higher-order terms over 1/N C ). 

As for galaxy surveys (Szapudi & Szalay 1997) we can generalize the estimators (13) to measure the cross-correlations 
between two surveys (or two subsamples) by defining: 



Mp «- (JVi)„(JV a ), 



iVi N 2 

/ , / J Qii e il ■ ■ ■ Qip£ipQji e jl ■ ■ ■ Qj q ( 

(ii,...,tp) O'li — iJa) 



(21) 



Obviously, these estimators could be extended to s— point correlations between s surveys. In this paper we shall only consider 
the case s i = s2 . Then, as in eq.(18) we can define the estimators M pq by averaging M pq over JV C cells. Of course, the 
one-point estimators can also be written M p = M p o or M p = Mo p depending on which survey they apply to. Next, we can 
define the analogs of the one-point estimators H p and H p . For the aperture mass, there is only one new cross-correlation at 
third order, M 2 \ (and its symmetric M12), and we define: 

H21 = Mn - 2MnM 10 - M20M01 whence (H21) = <M a 2 pl M ap2 ) c . (22) 

For the shear components 7i S we have two new cross-correlations at fourth order, M31 (and its symmetric M13) and M22, and 
we define: 

H 31 = M31 - 3.M11M20 -3X 20 Mii+ 3X11X20, whence (H31) = (7» 3 s i7is2>c = {llilis2) - 3(7iLi>(7«i7i82>, (23) 

and: 

#22 = M 22 - M 20 M 02 - M02M20 - 4MiiM n + M20M02 + 2X11X11, (24) 

whence: 

(H22) = {-yiiji^c = {^Lili^ - <7 l 2 si)<7 l 2 S 2) - 2(7 ls i7 is2 ) 2 . (25) 

Again, one can easily check that the scatter of H pq is always smaller than the scatter of M pq (see Appendix). Finally, from 
these estimators H p and H pq we can evaluate the parameters S p , t p and the cross-correlation coefficients r pq introduced in 
cqs.(ll), (12), through the estimators S p , T p and lZ pq : 

^^xf 1 ' % = ntv and Upq = h^/hI'%^' (26) 

* * ' L p+q,() ' L 0,p-\-q 
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We can note from numerical computations that the scatter of the moments (XjXf } increases very rapidly with the order 
p + q. Therefore, we can neglect the contributions to the scatter of S p due to the denominator. Hence we write: 



(S p ) ~ S p 

and: 
a 2 (TZ pq ) 



„(c \ ^ a C H p) 

& (Hpo) a (Hop) 
(X{)1 + (Xl)l 



v 2 (7-(-p, q ) 

(x{xi)l 



+ 



(T p ) ~ t p and 
a (TLpo ; Ti 



lap) 



(X{) C (X%) C 



p + q 



o\n 



v+qft) 



p+q\2 



{xr q ) 



(Kpi) 



(27) 
(28) 



p + q 



& CHQ,p+q) _ 2p a (Hp, q ;H p + q ,o) 

<Xf+ 9 } 2 p + g (XfAf > c (Xf +9 > c 



2g cr (Hp, q ;Ho,p+q) 2pq a (H p + q ,o;'Ho,p+ q ) 



p + <?(XfX|) c <Xf+«} c (p + g) 2 (Xf +9 ) c (Xf + «) c 



(29) 



Here we assumed that the scatter of T p and lZ pq is small so that it can be linearized in eqs.(26). If this is not the case, that is 
a2 (Tp) ^ 1 or a 2 (lZ pq ) <) 1, the error-bar obtained from eq.(28) or eq.(29) is not correct (it significantly underestimates the 
scatter) but this means that the measure is dominated by the noise so that we do not need an accurate estimate of a 2 : we 
only wish to be aware that the noise is too large to allow a precise measure of t p or r pq . In eqs.(28)-(29) we introduced the 
cross-dispersions a (Ti.p, q ;'Hp l , q l ) defined by: 

a (Hp, q ;H p \ q ') = (H P , q Hp\ q ') — {Hp, q ){Hpi, q i}. (30) 

We give in appendix A the expression of the various dispersions a 2 which we need in this article. 



4 NUMERICAL RESULTS 

In this section we numerically evaluate the signal to noise ratios involved with the various estimators described in the previous 
section. We assume a specific model for the background cosmology, a specific correlation hierarchy for the matter distribution 
and we focus on the SNAP observational strategy for numerical calculations. However as mentioned before the basic formalism 
developed here remains completely general and specific details studied here only serve illustrational purposes. 

4.1 Cosmological Parameters 

For the background cosmology we consider a fiducial LCDM model with Q m = 0.3, £Ia = 0.7, Ho — 70 km/s/Mpc and 
as — 0.88. To check the sensitivity of cross-correlations on the background dynamics of the universe we shall also study the 
effect of small variations of these parameters onto weak-lensing observables and their cross-correlations. We consider many- 
body correlations of the matter density field which obey the stellar model (Valageas et al. 2004; Munshi et al. 2004). This 
is actually identical to the minimal tree-model up to third-order moments (Munshi et al. 2004). This is also coupled to the 
fit to the non-linear power-spectrum P(k) of the dark matter density fluctuations given by Peacock & Dodds (1996). Lowest 
order non-Gaussian signatures will be at third order for M ap and at fourth order for smoothed shear components 7j s . 

4.2 Survey Parameters 

Hereafter, we adopt the characteristics of the SNAP mission as given in Refregier et al.(2004). More precisely, we consider 
the "Wide" survey where the redshift distribution of galaxies is given by: 



n(z s ) oc z s e (Zs ' zo) and zq = 1.13, 



The variance in shear due to intrinsic ellipticities and measurement errors is a* = (|e*| ) 



2\l/2 



(31) 

0.31. The survey covers an 

area A — 300 deg 2 and the surface density of usable galaxies is n g — 100 arcmin -2 . Therefore, we take for the number N of 
galaxies within a circular field of radius 9 S : 

n a 



N = n g Txei ~ 314 



lOOarcmin 2 J V larcmin 
and for the number N c of cells of radius 9. 



N c 



A 



(29s 



= 2.7 x 10 5 



A 



300deg' 



■)( 



larcmin 



(32) 



(33) 



For the shear this number somewhat overestimates JV C because of the sensitivity of *y iB to long wavelengths, which would 
require the centres of different cells to be separated by more than 29 3 in order to be uncorrelated. 
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Figure 2. Left panel shows the quantity p(0 B ) which measures the relative importance of galaxy intrinsic ellipticities, for the aperture 
mass and two redshift bins, as a function of the angular scale 9 S (in arcmin). The variance {M% ) for these two redshift bins is displayed 
in middle (low-z subsample) and right (high-z subsample) panels. Various line styles correspond to various cosmological or survey 
parameters. The solid "A" line is our fiducial model described in sections 4.1-4.2. The dotted "f2 m " curve shows the effect of a 10% 
increase of H m (from Q m = 0.3 up to f2 m = 0.33), the dot-dashed u crg" line corresponds to a 10% increase of ag (from ag = 0.88 up to 
ag = 0.97) while the dashed "zo" curve represents a 10% decrease of the characteristic redshift zq of the survey (from zq = 1.13 down 
to zq = 1.02). These line styles will be used in the following plots. The error bars show the one a dispersion around the mean value. 
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Figure 3. The quantities *2 and ru (see text for description) for Af ap are plotted as a function of smoothing angle 9 S for the two redshift 
bins. Error bars correspond again to the one cr dispersion around the mean value. 



In order to extract some information from the redshift dependence of weak lensing effects, which could also be used to 
discriminate noise sources, we also divide the "Wide" SNAP survey into two redshift bins. Then, we shall study the cross 
correlations between both bins. Thus, we consider the two subsamples which can be obtained from the "Wide" survey by 
dividing galaxies into two redshift bins: z s > z* and z s < z*. We choose z* = 1.23, which corresponds roughly to the separation 
provided by the SNAP filters and which splits the "Wide" SNAP survey into two samples with the same number of galaxies 



(hence 



50 arcmin ). Note that one cannot use too many redshift bins at it decreases the number of source galaxies 



associated with each subsample (for the aperture mass we could still obtain good results with three bins but we shall restrict 
ourselves to two redshift bins in this paper). The redshift bins that we use are similar to those which are used by Refregier 
et al.(2003) using photometric redshifts, except that they have a sharp cutoff and non-overlapping source distributions. Note 
that using overlapping source distributions (over redshift) would increase the cross-correlations. 



4.3 Numerical evaluation of Estimators and their scatter 

For all our computations we use the matter correlation hierarchy introduced and tested in a recent series of papers for the 
evaluation of statistics related to shear and M ap statistics. These analytical models were tested extensively against numerical 
simulations and were found to provide satisfactory results. We have studied the cross correlations among the two redshift bins 
obtained from the "Wide" SNAP survey using both the aperture mass M ap and the smoothed shear components 7^. In the 
following, in weak-lensing observables of the form H pq the first index p refers to the low-z subsample and the second index q 
to the high-z subsample. 

4-3. 1 Aperture mass statistics Af ap 
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Figure 4. The skewness 53 of the aperture mass is displayed as a function of the smoothing angle 9 S , for the low-z subsample (left 
panel) and the high- z subsample (right panel). 



4.3.1.1 Second-order cumulants We first present in the left panel of Fig. 2 the quantity p{6 s ) as a function of the 
angular smoothing radius 9 S for the two redshift bins obtained from the "Wide" SNAP survey. We can note that p is of 
order unity which means that the galaxy intrinsic ellipticities yield a significant contribution to the signal which cannot be 
neglected (but weak lensing effects can still be extracted) . The relative importance of intrinsic ellipticities decreases at large 
scales (p increases) as the number N of galaxies within the angular radius 9 B grows (this effect dominates over the decrease 
of the amplitude of weak lensing distortions, see eq.(16)). The quantity p is larger for the high-redshift subsample because its 
variance (Mf p ) c is larger (see eq.(16)). 

Next, we show in the middle and right panels of Fig. 2 the variance (Mf p ) of the aperture mass for each redshift subsample. 
The four curves correspond to our fiducial LCDM model and to a 10% variation of f2 m , as or zq. We can see that the error-bars 
are quite small and would allow an accurate measure of the cosmological parameters. However, the signal is also sensitive to 
the redshift distribution of the sources which may be an important limiting factor (although the dependence is rather weak). 
As is well known, the amplitude of weak lensing distortions increases with the redshift of the sources as the line of sight is 
longer, hence the variance is larger for the high-z subsample. On the other hand, it decreases at larger angular scales where 
the amplitude of density fluctuations is smaller. Note that the error-bars for the variance are quite small which justifies the 
approximation (27) for the scatter of S p . 

Then, we display in the left panel of Fig. 3 the ratio £2 between the variances (Aff p ) of the weak lensing effects associated 
with each redshift bin, as defined in eq.(ll). In agreement with Fig.2 we have £2 < 1. We can see that the sensitivity 
to cosmological parameters is rather modest. The error-bars are of the same order as the variation with a 10% change of 
cosmological parameters so that by combining several angular scales one can get a useful constraint on cosmology. In particular, 
note that fi m and as affect the ratio £2 in opposite directions while they affect the variances (Mf p ) in the same direction. Thus, 
dividing weak lensing surveys like the SNAP mission into several redshift bins provides additional information in addition to 
the overall amplitude of weak lensing effects averaged over all source redshifts. 

Finally, we present in the right panel of Fig. 3 the cross-correlation coefficient rn (defined in eq.(12)) between the two 
redshift bins, for our four models. We can see that the cross-correlation is quite strong (rn ~ 0.9) with a small scatter. 
Indeed, the lines of sight associated with the two subsamples probe the same density fluctuations at low z, where they are 
largest because of the build-up of large-scale structures. Moreover, the cross-correlation coefficient rn only shows a very weak 
dependence on cosmology or the redshift distribution of sources: as the two subsamples are highly correlated most of the 
dependence on cosmology or source redshifts which appeared in Fig.2 cancels out between the numerator and denominator in 
eq.(12). This also leads to reasonably small error-bars since the dispersions of the numerator and denominator do not merely 
add up in quadrature, see eq.(29). The fact that the dependence on cosmology is very weak (and smaller than error-bars) 
means that the cross-correlation coefficient rn and the cross-product (M a piM a pa) are not useful to constrain cosmological 
parameters. They do not add much information about cosmology to the one already included in the one-point statistics 
(M^pi) and (Mf p2 ). However, we can take advantage of this property to use rn to measure the observational noise. Indeed, 
since rn can be predicted up to a very good accuracy (e.g. from other statistics like the one-point variance within the same 
survey, or from other sources like the CMB) its value gives a simple estimate of the scatter of second-order cumulants as 
measured in a specific survey. It could also help track down systematics. As noticed in section 3, eq.(29) for the scatter of 
the cross-correlation coefficient rn assumed that the dispersion was small a -C 1. If this is not the case, that is the value rn 
obtained from a weak-lensing survey shows a large departure from its theoretical mean (e.g., one measures rn Ss 0.5), then 
eq.(29) is an underestimate. Besides, the probability distribution of rn would be strongly non-Gaussian. Then, the estimate 
of the dispersion a of second-order cumulants one would derive from the measure of rn would be too large. Therefore, the 
break-up of eq.(29) in this domain is not a real problem since it merely yields conservative estimates for the noise (besides in 
such cases the data is too noisy to provide accurate constraints on cosmology). 



© 0000 RAS, MNRAS 000, 000-000 



Cross-correlating Weak Lensing Surveys 9 




u 



8 10 



1 
0.95 

0.9 
0.85 

0.8 
0.75 




4 6 8 



~ 0.9 r 



u 




8 10 



Figure 5. The ratios t-$, r%\ and ri2, for the M ap statistic and two redshift bins, are plotted in left, middle and right panel respectively. 
Error bars correspond again to one a scatter. Note the very weak dependence on cosmology or source redshift. 

4.3.1.2 Third-order cumulants The variance of weak lensing observables cannot simultaneously constrain several cos- 
mological parameters. Thus, in order to break the fi m — as degeneracy it was proposed to consider the skewness S3 (Bernardeau 
et al. 1997). More generally, taking into account higher-order cumulants brings further information which can help constrain 
the underlying cosmology. However, the noise increases very rapidly with the order of the observed cumulants so that we shall 
restrict ourselves to the lowest-order non-vanishing cumulants beyond second-order. For the aperture mass this simply corre- 
sponds to third-order statistics like the skewness 5*3 defined in eq.(ll) (note that another possibility is to look at the whole 
probability distribution itself, see Valageas et al. 2004). We present in Fig. 4 the skewness S3 of the aperture mass as a function 
of the smoothing angle 9 S for the two redshift bins. In agreement with previous works we find that the skewness increases for 
lower source redshifts (because of the smaller length of the line of sight and of the growth with time of the non-Gaussianities 
of the density field). Hence we can check that 530 > 5o3- We can see that although the scatter is much larger than for the 
variance, one should obtain a clear detection of non-Gaussianity in all cases. Moreover, by combining different angular scales 
and redshift bins (e.g. through a Fisher matrix approach) one should be able to constrain cosmological parameters at a 10% 
level (note however the large dependence on the source redshift, which is more important than for second-order cumulants). 

We present in the left panel of Fig. 5 the ratio £3 between the third-order weak-lensing cumulants associated with the two 
redshift bins. In agreement with Fig. 4 we obtain £3 < 1. We note that the dependence on cosmology is very small and much 
below error-bars. Therefore, redshift binning is not very useful to constrain cosmological parameters when it is applied to 
these third-order statistics. On the other hand, the measure of £3 could be used to evaluate the noise of third-order cumulants 
(or the skewness) since its value can be predicted with a good accuracy (without requiring a high accuracy for cosmological 
parameters). This can be useful to assess the reliability of the data. 

Next, we display in the middle and right panels of Fig.5 the cross-correlation coefficients r^\ and rvx between both 
subsamples, as defined in eq.(12). We can check that the cross-correlation is very strong (r ~ 0.9) which means that the weak 
lensing effects associated with the two redshift bins are highly correlated. Of course, this is consistent with the high value 
already obtained for the cross-correlation m of second-order cumulants (Fig. 3). Note that we get r^\ > ri2 > r\\. Indeed, 
higher-order cumulants give more weight to the common row-z parts of the line of sights (because the parameters S p of the 
3-d density field increase in the non- linear regime). This increases the correlation of higher-order weak lensing cumulants. 
Similarly, the coefficient rzi is slightly larger than m as it gives more weight to low- 2 distortions. We can note that the scatter 
of the cross-correlation coefficients r<z\ and ri2 is somewhat smaller than might be guessed from the dispersion of third-order 
cumulants (Fig. 4). This is again due to their high correlation. As for the second-order cross-correlation rii, the dependence 
on cosmology or the source redshifts is very small (and much below error-bars). This means that r2i and ri2, as the cross- 
products (Mf pl M ap 2)c and (M ap iAf a ! p 2}c, do not bring additional information to the one-point cumulants S3 with regard to 
cosmology. However, even though they are mostly useless with respect to the measure of cosmological parameters, they can 
be very useful to evaluate the noise of third-order cumulants. Indeed, a simple measure of r2i or 7-12 allows a good estimate 
of the noise of the skewness since their theoretical value can be predicted with a good accuracy. Therefore, cross-correlations 
are a complementary tool to one-point statistics and should prove useful. 

4-3.2 Smoothed shear components 7i s 

We now describe our results for the shear components, following the outline of our study of the aperture mass. 

4.3.2.1 Second-order cumulants We first show in the left panel of Fig. 6 the quantity p which measures the importance 
of galaxy intrinsic ellipticities, while we show the variance (7^) in the middle panel (low-z bin) and the right panel (high-z 
bin). As is well known, the variance of the shear components is larger than the variance of the aperture mass because the 
compensated filter of M ap damps the contributions from long wavelengths. This also leads to a larger value of p for the shear 
components. Hence the relative size of the error bars associated with second-order statistics are smaller for the shear. The 
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Figure 6. Same as Fig. 2 but for shear components. Left panel shows the characteristic function p associated with both redshift bins. 
Middle and right panels show the variance (7? ) as a function of smoothing angle 8 S for lower and higher redshift bins respectively. 
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Figure 7. Same as Fig. 3 but for shear components. As in case of Af ap , £2 and rii are almost independent of cosmological parameters. 

magnitude of these effects depends on the slope of the matter density power-spectrum and they are larger at smaller angular 
scales. Again, we can check that the error bars are very small and the variance associated with each bin should be accurately 
measured. This would yield a strong constraint on cosmology. 

Next, we present in the left panel of Fig. 7 the ratio £2 of both variances (7^) associated with the two redshift subsamples. 
Although the dependence on cosmology is rather weak, it could be extracted from the data. Contrary to the aperture mass, 
fljn and erg affect the ratio £2 in the same direction as each variance (7? s } but with different powers. Therefore, redshift binning 
still provides additional constraints on cosmology. Finally, we display in the right panel of Fig. 7 the cross-correlation coefficient 
m. As for M ap it is rather high, which translates a strong correlation between both weak-lensing signals, and it only shows 
a very weak dependence on cosmology. However, since error bars are quite small it might be possible to use rn to further 
constrain cosmological parameters. On the other hand, this also means that it does not provide a very robust tool to estimate 
the noise of the survey, unless the cosmological parameters and the source redshifts are known up to a high accuracy. 

4.3.2.2 Fourth-order cumulants As for M ap , we now study the lowest-order cumulants which probe the deviations 
from Gaussianity. Since the shear components are even random variables we need to go up to fourth order. Since the noise 
increases very fast with the order this means that non-Gaussianities are more difficult to detect from the shear components 
than from M ap (where third-order cumulants do not vanish). Thus, we display in Fig. 8 the kurtosis S4 of the smoothed shear 
components for both redshift subsamples. As for the skewness of M ap , the kurtosis increases for lower source redshifts, and we 
obtain S40 > So4- As expected, we can check that the error bars are much larger than for the skewness of M ap . In particular, 
at large angular scales (beyond 6') it will be difficult to obtain a clear detection of non-Gaussianity. At lower angular scales 
one should be able to obtain a meaningful measure of S4. On the other hand, one should note that the theoretical error bars 
for S4 are actually rather large in this transition regime towards non-linear gravitational structures, so that one cannot expect 
to obtain strong constraints on cosmology from 1S4. Rather, it is probably safer to use S4 to check the broad consistency of 
the gravitational clustering process. 

Next, we show in the left panel of Fig. 9 the ratio £4 between the fourth-order cumulants associated with both redshift 
subsamples. As for M ap , we get £4 < 1 and the dependence on cosmology is very small. Again, this means that redshift binning 
does not provide much additional information for fourth-order statistics. On the other hand, a measure of £4 could provide a 
convenient estimate of the noise of the survey with regard to these fourth-order cumulants. 

Finally, we present in the middle and right panels of Fig. 9 the cross-correlation coefficients r3i,ri3 and r22- We obtain 
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Figure 8. The kurtosis of smoothed shear components S4 is plotted for two redshift bins as a function of smoothing angle 9 S (low-z 
subsample: left panel, high-z subsample: right panel). 




2 4 6 8 10 



S- 



1.1 



0.9 



0.8 



2 4 6 8 10 
0- 



f-t 



1 . 1 










1 


It 










- : 


0.9 












— 


0.8 












-~ 




-i i i i I i i 


I 


I I 


i i - 



2 4 6 8 



Figure 9. The left panel shows £4 as a function of smoothing angle 9 S . The middle and right panels correspond to r3i,ri3 and r22 
respectively. These quantities are related to the kurtosis of the shear field and are dominated by the noise. We have reduced the error 
bars by a factor of four for plotting. 

rm > ri3 > r-zz- Again, the cross-correlation is large and the dependence on cosmology is very weak and much smaller than 
error bars. Therefore, these cross-correlations are mostly useful to estimate the noise of the survey rather than the cosmological 
parameters. We must note that at large angular scales the formula (29) for a 2 (1Z pq ) yields a negative number (in which case 



with plot a 



2,1/2 



This translates the fact that error bars are so large that errors cannot be linearized in the expression 



(26) of IZpq. However, in these cases we have a <; 1 so that an accurate estimate of the scatter is not needed: the signal is 
completely dominated by the noise and useless for practical purposes. 

5 DISCUSSION AND FUTURE PROSPECTS 

In order to derive precise results from weak lensing surveys one must take into account various sources of noise like the intrinsic 
ellipticity distribution of galaxies and the cosmic variance. Extending a previous study (Valageas, Munshi & Barber 2004) 
we have introduced these realistic sources of noise in various estimates of cross-correlation statistics between several redshift 
subsamples. Such studies should help discriminating weak lensing effects from noise and measuring cosmological parameters. 

Thus, extending an earlier work (Valageas, Munshi & Barber 2004) where we focused on one-point statistics, we have 
taken advantage of the high projected density of galaxies n g in the future SNAP survey to divide the source population into 
two redshift bins. Then, we have computed cross-correlations between these subsamples. Note that our formalism is kept 
completely general and can be applied to different surveys with or without any overlap in source redshift. We have used two 
statistics, the aperture mass M ap and the shear components ji s , to quantify the correlations among redshift bins as a function 
of the smoothing angle 9 S . The underlying model that we have used for computing the cumulant correlators is the same as 
that of Valageas, Munshi & Barber (2004). This has been tested against numerical simulations for a wide range of smoothing 
angles and source redshifts and was found to provide good results. 

We have considered both second-order cumulants (cross-correlation rn) and the lowest-order cumulants which can mea- 
sure deviations from Gaussianity (third or fourth order for M ap or 7j b ). High order cumulants can remove the well-known 
degeneracy between S2 m and as but the relative importance of the noise increases rapidly with order. This makes the aperture 
mass a better tool than the shear components which require fourth-order statistics. However, note that weak lensing statistics 
are also sensitive to source redshifts which might prove a significant limitation for any weak lensing survey. We noticed that 
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for a survey such as SNAP one could still obtain interesting information from the shear components after dividing the sample 
into two redshift bins, while for the aperture mass one could go up to three subsamples. Increasing further the number of 
subsamples increases the noise of third-order or higher order cumulants and is mostly useless (except for the variance where 
the error bars are small). We have found that cross-correlations between subsamples do not bring additional information with 
regard to the constraints on cosmology to those already provided by one-point statistics (except marginally for the second- 
order moments of the shear). Indeed, different source populations selected by their redshift (in the zero angular separation 
case) still probe the same density fluctuations at low z (where they are largest). This implies a high value for the various 
cross-correlation coefficients which are almost independent of cosmological parameters. These properties hold even better 
for higher orders statistics. However, this feature can be used to estimate the noise of the survey. Therefore, one-point and 
two-point statistics provide complementary tools: the former can constrain cosmology and check the gravitational clustering 
process, while the latter can yield an estimate of the noise of the survey. 

Recent studies have underlined the importance of joint estimation of cosmological parameters such as galaxy-shear cross 
correlation, shear-shear correlations along with galaxy-galaxy angular correlations. Deep multi-colour galaxy surveys with 
photometric redshifts will provide an unique opportunity to study such large numbers of two-point correlation observables to 
pin-point the background dynamics of the universe as well as non-Gaussianity and bias associated with cosmological density 
distribution of underlying mass as well as galaxy populations. The statistical machinery developed here will be extremely 
useful in this direction. 

Although the cross correlations and the error terms are defined in this work with a specific filter in mind and we only 
concentrate on two different statistics the framework developed here has been kept completely general and the formalism can 
be used to compute the signal to noise related to quantities such as foreground-background galaxy populations (Moessner & 
Jain 1998). Fluctuations in size distribution of source galaxies can also be used to study the weak lensing convergence directly 
without having to go through the map making procedure, from shear maps to convergence maps (see e.g. Jain, 2002). Higher 
order cross-correlation statistics or cumulant correlators we develop here can be useful in probing correlations of such measures 
of weak lensing convergence with cosmic shear measurements using shapes of background galaxies as observations improve 
and surveys start to cover a larger part of the sky and probe deeper. Not only these statistics can provide valuable consistency 
checks among various measures from different surveys, they will also be a useful supplement in constraining cosmological 
parameters when used along with their one-point counterparts. 

Recently Dalai et al. (2002) have investigated the correlation between the smoothed convergence and the value along 
the line of sight to individual supernovae. However these calculations were done by ignoring non-Gaussian corrections. The 
formalism developed here takes into account all non-linear corrections from higher order contributions. In our formalism we 
have shown that such calculations can be treated as a special case of the generic expression derived here in a more general 
context. These results can not only be used to analyze cross-correlations among various redshift surveys but also to probe 
supernova weak-lensing cross-correlation. A more detailed analysis will be presented elsewhere. 

Although current studies of weak lensing effects have focused on the distortion of background galaxy ellipticities, it has 
been pointed out recently that such studies can be augmented in the near future by lensing studies of unresolved sources 
and also of the spatial fluctuations of their integrated diffuse emission. In particular one could study the diffuse background 
generated at far-infrared wavelengths by dusty star bursts, first stars and galaxies in near-infrared and also the 21cm emission 
from neutral gas which forms the intergalactic medium prior to reionization. These studies promise to smoothly interpolate 
from the weak-lensing studies of the cosmic microwave background at very high redshift downto weak lensing studies using 
galaxy shapes at low redshifts. Cross-correlation studies of one type of surveys with another hold the prospect of mapping the 
dark matter distribution of the universe at medium redshift. We plan to present results of such analysis in a separate paper. 

In almost all weak lensing studies one uses the Born-approximation. Its validity is checked in numerous numerical studies 
of weak lensing at small angular scales. In addition we have ignored contributions from source clustering and lens coupling. 
Some of these issues have been studied in Bernardeau (1997), Bernardeau et al (1997) and Schneider et al. (1998). It was 
found that at large smoothing angular scales these contributions are negligible. In the highly nonlinear regime an accurate 
picture of galaxy bias is required to deal with these issues but it is still lacking at the moment. 

The correlations of galaxy intrinsic ellipticities might yield additional complexities in dealing with shear correlations and 
elaborate schemes have been developed (Heymans & Heavens 2003, Crittenden et al. 2002). However the extent to which such 
correlations will affect weak lensing surveys remains somewhat uncertain. It is generally believed that such effects will play a 
less important role as we increase the survey depth and we can reduce their role through acquisition of photometric redshift 
(Heymans & Heavens 2003). 

It is almost always assumed in all weak lensing studies that the galaxy intrinsic ellipticity distribution is Gaussian although 
it might not be the case. Any signature of such non-Gaussianity if found by observational teams will have to be included into 
our analytical calculations. However this can be performed in a straightforward way in our formalism. Clearly, if such intrinsic 
non-Gaussianities are too large they might even dominate the signal and complicate inferring the observational data. 

It became clear that future weak lensing surveys will mostly be limited by sky coverage which determines the scatter due 
to the finite size of the catalog. The intrinsic ellipticities will only dominate at small angular scales and will not preclude a 
meaningful estimation even with small sky coverage. In our analysis we have used the plane parallel approximation to compute 
the higher order non-Gaussianities. However as the survey size increases in the future it will be important to translate these 
studies using all-sky calculations. A detailed analysis will be presented elsewhere. 
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APPENDIX A: SCATTER OF LOW-ORDER ESTIMATORS M PQ AND H PQ 

We give in this appendix the expression of the scatter a 2 of the various estimators M vq and H pq introduced in section 3 which 
we need to derive the numerical results presented in section 4. Note that although in our numerical calculations we take the 
smoothing windows 9 S of different surveys or subsamples to be the same when we compute cross-correlations, this needs not 
to be the case for the following equations to be valid. 
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Al Aperture Mass 

We first consider the aperture mass statistics, where we restrict ourselves to cumulants of order two and three. The dispersions 
a of first and second-order estimators are (with (Mi) = 0): 



1 + i 
P 



a 2 (M 2 ) = (M 4 p ) c + (M a 2 p } 2 2 



a\M,) = (M4) c 
and: 

a 2 (M n ) = <M a 2 pl M 2 p2 ) c + (M apl M ap2 > 2 + (M a 2 pl ) c (M a 2 p2 ) c 



1 + i 
P 



1+1" 

P l . 




1+1 

pi 



(Al) 



(A2) 



Note that in the case p\ = pi the cross-correlation error <j 2 (Mh) does not become identical to the error <j 2 (M 2 ) associated 
with one survey for the estimator M 2 . This simply reflects our assumption that the intrinsic ellipticities of different galaxies 
are not correlated and that different surveys or subsamples do not share any galaxy. Of course, it is also possible to handle 
the case where a certain fraction of the source population is shared by both surveys but we shall not consider this case here. 
As seen in eqs.(29), (30), we also need the dispersions of cross-products: 



a 2 (Mn;M 2 o) = (M a 3 pl M ap2 ) c + (M a 2 pl ) c (M ap iM ap2 ) c 2 



1 + 1 
Pi 



a 2 (M 20 ; M 02 ) = <M a 2 pl M 2 p2 ) c + 2(M apl M ap2 ) 2 . (A3) 



Note that the galaxy intrinsic ellipticities do not contribute to the dispersion of the product (M20M02) because we assumed 
that the two subsamples have no common galaxies. Next, the scatters of third-order estimators are: 



a 2 (M 3 ) = (M a b p ) c + (M a 4 p ) c (M 2 p ) c 
and: 



15 + ? 
P 



+ 9<M 3 p > 2 + (M a 2 p ) 3 



lr 27 18 6 

15 H h-r + ^r 

p p z p A 



(A4) 



a(M 21 ) = <M a 4 pl M a 2 p2 ) c + <M a 4 pl ) c (M aP . 2 ) c 



1 + 1 
Pi 



+ 8(M a 3 pl M ap2 ) c (M ap iM ap2 ) c + (M a 2 pl M 2 p2 ) c (M a 2 pl ) c 



6 + 1 
Pi 



+4(M a 3 pl ) c (M apl M a 2 p2 ) c + 5(M a 2 pl M ap2 ) 2 + (M a 2 pl ) 2 (M a 2 p2 ) c 

8 



n 4 3 

3 + — + — 



+ 



Pi P'2 P'{ PlPl P 2 iP'2 



+(M ap iM ap2 ) 2 (M a 2 pl ) c 



12 + 



Pi 



(A5) 



For cross-products we obtain: 



a 2 (M 2i ; M30) = (M a 5 pl M ap2 ) c + 5(M 4 pl ) c (M ap iM ap2 ) c + (M 3 pl M ap2 ) c (M 2 pl ) c 



10 + 1 
Pi 



+ 9(M a 2 pl M ap2 ) c (M a 3 pl ) c 



+(M ap iM ap2 ) c (M a 2 pl ) 2 



lr 18 6 

15 + — + ^T 
Pi Pi 



(A6) 



a 2 (M 21 ; M 0a ) = (M 2 pl M 4 p2 ) c + (M 4 p2 ) c (M 2 pl ) c + 8(M ap iM 3 p2 ) c (M ap iM ap2 ) c + (M 2 pl M 2 p2 ) c (M 2 p2 } c 



6 + 1 

P'2 



+3(M a 2 pl M ap2 ) c (M a 3 p2 ) c + 6(M apl M a 2 p2 ) 2 + (M a 2 pl ) c (M a 2 p2 ) 2 



3+1 

P2 



+ (M apl M ap2 ) 2 (M 2 p2 ) c 



12 + 1 
pi 



and: 



(A7) 



(A8) 



a 2 (M 30 ; M 03 ) = <M a 3 pl M a 3 p2 ) c + 3(M a 3 pl M ap2 ) c (M a 2 p2 ) c + 3(M ap iMf p2 ) c (M a 2 pl ) c + 9(M a 2 pl M a 2 p2 ) c (M apl M ap2 ) c 
+9(M 2 pl M ap2 ) c (M apl M a 2 p2 ) c + 9(M a 2 pl ) c (M a 2 p2 ) c (M apl M ap2 ) c + 6(M ap iM ap2 ) 3 . 

As recalled in section 3, following Valageas et al. (2004), it is better to use the cumulant estimators H pq rather than the 
moment estimators M vq , as they show a smaller scatter. Their dispersion can be expressed in terms of the scatter of the 
estimators M pq as: 



a 2 (H 3 ) = a 2 (M 3 ) - 6(M 4 p ) c (M 2 p ) c - 9(M a 2 p ) 3 



1 + i 
P 



(A9) 



a 2 (H 21 ) = a 2 (M 21 ) - 4(M a 3 pl M ap2 ) c (M apl M ap2 ) c - 2(M a 2 pl M a 2 p2 ) c (M a 2 pl ) c - (M 2 pl ) 2 (M a 2 p2 ) c 



1 + 1 
Pi 



-(M apl M ap2 ) 2 (M a 2 pl ) c 



8 + 1 
Pi 



(A10) 
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« 2 (H 21 ;H 30 ) = a 2 (M 21 ; M 30 ) - 2(M a 4 pl ) c (M apl M ap2 ) c - 4(M a 3 pl M ap2 ) c (M a 2 pl ) c - (M apl M ap2 ) c (M a 2 pl ) 2 



9 + 



Pi 



(All) 



a 2 (H 21 ;H 03 ) = a 2 (M 21 ; M 03 ) - <M 4 p2 ) c <M a 2 pl ) c - 2(M apl M a 3 p2 ) c (M apl M ap2 ) c - 3<M a 2 pl M 2 p2 ) c (M a 2 p2 ) c 



-<M a 2 pl ) c <M a 2 p2 } 2 



3 + A 

P2 



-6(M apl M ap2 ) 2 (M a 2 p2 ) c 



(A12) 



and: 



a 2 (H 30 ;H 03 ) = a 2 (M 30 ; M 03 ) - 3<M a 3 pl M ap2 ) c (M a 2 p2 ) c - 3<M apl M a 3 p2 ) c (M a 2 pl ) c - 9(M a 2 pl ) c (M a 2 p2 ) c (M apl M ap2 ) c . (A13) 

We can check in eqs.(A9)-(A13) that the scatter of the estimators H pq is smaller than the dispersion of the estimators M vq . 
Besides, one can note that the galaxy intrinsic ellipticity dispersion only enters the expressions of the scatter of the estimators 
H pq through the combination (1 + 1/p). 



A2 Shear components 

The dispersions of second-order cumulants are given by the same expressions (A1)-(A3) as the one obtained for the aperture 
mass. On the other hand, since all odd-order moments of shear components vanish we need to go up to fourth-order moments 
in order to obtain the first measure of deviations from Gaussianity. Thus, we obtain for fourth-order statistics: 



a 2 (M 4 ) = <7? s )c + <7? s )c<7&>c 



28 + ^ 
P 



+ 34<7 4 s ) 2 + < 7 4 } c (7 2 } 2 



nrtA 240 72 

204 + + — 

P P 2 



+ {lis)c 



„„ 240 216 96 24 

96 + + — + — + — 

p p 2 - p A p 4 



(A14) 



0" 2 (M 3 i) = (7fsl7» 2 s2)c + (7fsl>c(7?92>c 



1 + 



pi 



+ 12(7f s i7is2)c(7isl7is2)c + {'yil^} c{~fi B i) c 



15- 



pl 



+15(7. 4 si)c(7 l 2 3i7 l 2 32>c + 19(7 3 sl 7i s2 } 2 + (7fci)c(7?si)c(7?s2>c 



9 15 9 

15 + — + — + 



+ (7 l 3 al7«s2) C <7isl7is2} C (7 I 2 sl) C 



72 

114 + — 
Pi 



2 2 \ / 2 \2 



+ <7isl7is2>c(7i s l) 



pi p2 plp2 



54 18 



+ 30(7j 4 s i} c <7isi7»s2) 2 



45 H h — 

Pi Pi 



+ (7isl)c(7is2)c 



1K 27 15 18 27 6 18 6 

15 H 1 h^H h -» + —, 1- -s— 

Pi p2 p{ pip'2 pf p(p2 p{p2 



+ (7»sl7is2) c (7isl)c 



81 



108 36 
pi p{ 



,(A15) 



and: 



6 + ± 
Pi 



+ (7fel7. ? s2)c(7iU)c 



cr 2 (M 22 ) = (7 4 1 7 4 2 } c + 16(7 3 sl 7 3 s2 } c (7 ls i7 ls2 ) c + {-ytiita) oilti) 

+ 16(7i S i7 3 2 ) c (7»s 2 7i 3 si)c + 17(7?si7iU)c + 68(7 2 s i7 2 s2 ) c (7 is i7»s2) 2 + (-yti) c{lis2)l 

4 2 

+ (7 I 4 a2)c(7 I 2 1} 2 3 H h — + (7 3 1 7i s2 ) c (7isi7»s2) c (7i ! S 2} 

[ Pi Pi J L 

+ <7 l 2 s i7 l 2 S 2} c <7 ! 2 si) c (7 l 2 s2}c 34 + — + — + — — + 20(7 ls i7 ls2 } 4 + (7 2 si}c(7 2 s2} c <7«si7«s2) 2 
+ (7« 2 s i)?<7 2 32} 2 



6+± 

P2 



+ (7 l 4 sl)c<7i 1 s2}c 



48+^ 

P2 



n 4 2 

3 + — + ^ 

P2 P 2 _ 

+ <7isi7i i s 2 ) c <7^i7^2} c <7 2 si} c 



32 
Pi 



pi P2 plp2 



„„ 48 48 32 

68 + — + — + 



Pl p2 plp2 



„ 12 12 6 6 16 8 8 4 

8H 1 ^ + ^H \-— 2 \-— 2 ^^-T 

P2 Pl Pi P 2 P1P2 pfp2 p 2 pl PiP 2 

Next, for cross-products we obtain: 



(A16) 



a (M 3 i;M 40 ) = (7isi7;s 2 ) c + 7(7 isl ) c (7 is i7; s2 ) c + (7isi7»s2)c(7isi)c 



21 + 



12 
Pi 



+ <7 3 si7*s2)c(7 2 s i} 2 

+ <7isl7is2)c(7 l 2 l)c 



lr> „ 120 36 

102 + + -j 

Pi Pi 



+ (7isl)c<7isl7is2}c(7isl} 



+ 34(7 3 1 7 is2 ) c (7 4 1 ) c 
60 " 



102 + 



Pi 



96- 



180 108 24 
Pi P? P? 



(A17) 



CT (M 3 i; M 04 ) = <7i S l7» 5 s2)c + 3(7isl7» 5 s2}c<7isl}c + (7isl7is2)c<7 lS 2}c 
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+4(7i s i7; S 2} c (7is2}c + 30(7i s i7i s2 } c (7isi7i S 2) c + (7isi7is2)c(7*s2) 

12" 
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(A18) 



a (M40; M04) = (7isi7»s2)c + 6(7 isl 7 is2 } c (7 ls i) c + 6(7; s i7 is2 ) c (7 ls2 } c + 16(7i sl 7- s2 ) c (7 ls i7i S 2)c + 16(7 is i7is2}c(7i S i7 lS 2}c 

+18(7 J 2 a i7 l 2 S 2} 2 + 48(7isi7; , s2 ) c (7 2 si) c <7»si7is2} c + 48(7 3 s i7 ls2 ) c (7 2 s2 ) c (7 ls i7 ls2 } c + 36(7 2 s i7 2 2 ) c (7 2 s i) c (7 2 s2 } c 
+72( 7 2 sl 7 2 2 ) c (7 lsl7ls2 ) 2 + 72( 7lsl 7 ls2 ) 2 ( 7 2 sl ) c (7 2 2 ) c + 24(7; s i7 is2 ) 4 , (A19) 

and: 

8 



a 2 (M 22 ;M 4 o) = (7 l 6 s i7« 2 s2)c + <7fsi)c(7 2 S 2) c + 12(7!i7«<2)c(7isi7is2)c + {jii^) ci^isi) c 

+20(7, ? s i7i B 2>c + <7 l 2 i7?s2)c<7 2 si} 2 42 H h — + (lii) c(li B i) 0(7^2} 

L Pi Pi J 
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L P 1 } 

Finally, we obtain for the estimators H pq : 
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o- 2 (#2 2 ) = cr 2 (Af 22 ) - 8(7f sl 7 l 3 s2 ) c (7isi7is2)c - 2(7 2 1 7 4 2 ) c (7 2 1 ) c - 2(7* sl 7? 2 )c(7iU)c - 48(7 2 1 7 2 s2 ) c (7 ls i7i S 2^ 
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a (H 3 i;H4o) = o (M 3 i;M 40 ) -3(7 isl } c (7 ls i7i s2 ) c -9(7 l J sl 7i S 2} c <7 is i} c - <7 is i7i S 2} c (7isi)c 
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(A24) 



a (H 31 ;H 4) = o- (M 3 i;M 04 ) - 3(7 ls i7i 2 ) c (7 isl ) c - 6(7; sl 7 is2 ) c (7 is2 ) c - 3(7 lsl 7 is2 ) c (7; s i7*s2} c - <7i s i7;s2} c (7 ls i) c 
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(A27) 



As for the aperture mass, we can check that the scatter of these estimators H pq is smaller than for their counterparts M p 
and that the galaxy intrinsic ellipticities only enter their dispersion through the combination (1 + 1/p). 
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